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ABSTRACT 

The gravitational waves generated during supermassive black hole (SMBH) coales¬ 
cence are prime candidates for detection by the satellite LISA. We use the extended 
Press-Schechter formalism combined with empirically motivated estimates for the 
SMBH-dark matter halo mass relation and SMBH occupation fraction to estimate 
the maximum coalescence rate for major SMBH mergers. Assuming efficient binary 
coalescence, and guided by the lowest nuclear black hole mass inferred in local galac¬ 
tic bulges and nearby low-luminosity active galactic nuclei (~ 10® M©) we predict 
approximately 15 detections per year at a signal to noise greater than five, in each of 
the inspiral and ringdown phases. Rare coalescences between SMBHs having masses 
in excess of 10^ M 0 will be more readily detected via gravitational waves from the 
ringdown phase. 
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1 INTRODUCTION 

There is evidence for the existence of supermassive black 
holes (SMBHs) over a range of epochs; SMBHs with masses 
10® — 10® Mq are ubiquitous in local galactic bulges (e.g. 
Magorrian et al. 1998: Ferrarese 2002), while the SMBHs 
powering quasars [which have been detected at redshifts as 
high as 2 ; ~ 6 (Fan et al. 2001; Fan et al. 2003)] are es¬ 
timated to range between 10® — 10® Mq (see e.g. Shields 
et al. 2003). Evidence for the early assembly of SMBHs, 
when interpreted within a hierarchical cosmology, suggests 
that SMBH coalescence may be a frequent event. If this is 
the case, the gravitational waves (GWs) generated during 
SMBH coalescence are a prime candidate for detection by 
the Laser Interferometer Space Antennae (LISA^, Folkner 
et al. 1998; Flanagan & Hughes 1998; Hughes et al. 2001). 

Estimates of the SMBH coalescence rate depend cru¬ 
cially on the occupation fraction of SMBHs in halos, and 
therefore on the adopted model for the formation and growth 
of SMBHs. Within a hierarchical CDM cosmology, both 
seeded SMBH formation models [in which a primordial pop¬ 
ulation of SMBHs grow via accretion and/or merging (eg. 
Volonteri et al. 2003)] and ongoing formation models [where 
SMBHs form jointly with bulges in halos (eg. Kauffmann & 
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Haehnelt 2000)[ are consistent with the present-day ubiquity 
of SMBHs in galactic bulges (Menou et al. 2001). 

The formation of a SMBH may be limited to poten¬ 
tial wells above a minimum depth. If so, then there exists 
a critical halo virial temperature (Tvir.min) below which a 
halo cannot host SMBH formation. If Tvir.min is low (eg. 
Tvir.min ~ 10® K), then SMBHs are abundant in small 
halos, and predicted event rates are in the order of lOO’s 
per year (Wyithe & Loeb 2003; Haehnelt 2004). This event 
rate may be significantly lower, less than one per year, if 
Tvir.min ~ 10® K and SMBHs form only at the centres of 
massive galaxies (Haehnelt 2004). Hence the detection rate 
has the potential to constrain the global SMBH population 
[with the important caveat of the unknown efficiency of bi¬ 
nary black hole (BBH) coalescence, see jEl . 

The rate of coalescence detectable by LISA depends on 
the form of the GW signal and the instrument’s sensitivity 
curve. Previous estimates of the event rate have used charac¬ 
teristic (rather than time-dependent) GW signals to deter¬ 
mine approximate detection limits (Wyithe & Loeb 2003) or 
have estimated the number of detectable events by compar¬ 
ing the sensitivity curve with the total gravitational wave 
background due to coalescing BBHs (Sesana et al. 2004a). 
Recently Sesana et al. (2004b) estimated the expected event 
rate for a detailed physical model for SMBH growth. We 
perform the first empirically motivated calculations for the 
LISA detection rate of SMBH mergers in a hierarchical cos¬ 
mology that use accurate signal to noise ratios (SNRs) to 
determine detection criteria. In ^we describe the calcula- 
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tion of SNRs for the LISA detection of BBH coalescence. 
In ^we discuss our calculation of the SMBH merger rate, 
including our halo merger rate predictions, estimate for the 
occupation fraction of SMBHs in halos and the Mbh-Afhaio 
relationship. Finally, we present our event rate predictions as 
a function of Tvir.min O before concluding in ^ Through¬ 
out this work we assume flm = 0.27, Q,k ~ 0, Ha = 0.73, 
Ho = lOOh km s“^ Mpc“^ = 71 km s“^ Mpc“^, os = 0.84 
and a primordial power spectrum with slope n = 1 as 
determined by the Wilkinson Microwave Anisotropy Probe 
(WMAP, Spergel et al. 2003). 


2 BINARY BLACK HOLE COALESCENCE 

AND DETECTION BY LISA 

The efficiency with which BBHs coalesce is highly uncer¬ 
tain. Initially, the SMBHs sink independently toward the 
centre of a merged system due to dynamical friction from 
the dark matter background until they form a bound binary 
(Begelman, Blandford & Rees 1980). The efficiency of this 
process depends on the orbital parameters of the merging 
dark matter halos (van den Bosch et al. 1999; Colpi, Mayer 
& Governato 1999). As the orbital separation (r) decreases, 
3-body interactions with stars that pass within « r of the 
BBH centre of mass (the “loss cone”) increasingly dominate 
the energy loss. Depending on the orbital parameters of the 
binary and the background distribution of stars, this process 
may result in a hardened [orbital velocity (v) ~ stellar ve¬ 
locity dispersion (cr)] binary system. If hardening continues 
until a binary separation where energy losses are dominated 
by GWs, the binary will coalesce. 

Yu (2002) argues that the efficiency of BBH coalescence 
hinges on the time-scale for BBH hardening during the hard 
binary phase. During this stage, three-body interactions be¬ 
tween the BBH and individual stars eject stars into highly el¬ 
liptical orbits, lowering the inner stellar density and slowing 
further hardening (eg. Volonteri et al. 2003). The decelera¬ 
tion of the hardening rate is compounded by the preferential 
depletion of stars within the loss cone (region of parameter 
space where angular momentum is low enough that stars 
pass near the BBH), resulting in a deficiency of stars that 
can extract energy from the binary system (see Yu 2002). 

Uncertainties in the efficiency of processes which may 
replenish loss cone stars make this effect difficult to anal¬ 
yse. Two-body stellar relaxation is expected to result in 
some diffusion of stars back into the loss cone (Binney & 
Tremaine 1987) but it is unclear that this process alone can 
support sufficient hardening. Numerical N-body simulations 
have been used to include the effect of the BBH ’wander’ 
within the star held (due to 3-body interactions with stars), 
which increases the effective size of the loss cone and may 
prevent the coalescence from stalling (Quinlan & Hernquist 
1997; Milosavljevic & Merritt 2001; Chatterjee, Hernquist & 
Loeb 2003). Other scenarios which may aid evolution of the 
BBH into the GW dominated regime have been proposed, 
including the effects of gaseous disks (Gould & Rix 2000), 
battened or triaxial stellar distributions (Yu 2002) and dis¬ 
ruption by a third SMBH (Hut & Rees 1992). 

The efficiency of BBH coalescence is of prime concern 
for LISA detection rates. Whilst undoubtably still an open 
question, we note that there are plausible mechanisms to 


extract the required energy from the binary system. For the 
purposes of our calculation we therefore assume that 
hardening of the binary to the GW dominated regime occurs 
with an efficiency (emrg) of unity. All event rates quoted in 
this paper are proportional to emrg. 

The gravitational wave dominated portion of BBH 
coalescence may be divided into three main phases 
(Hughes 2002). The binary begins in the inspiral phase, 
with the SMBHs slowly spiralling into tighter orbits due 
to the adiabatic loss of GW energy. Eventually the dynam¬ 
ics become relativistically unstable and the SMBHs violently 
plunge to form a single object [merger phase). The hnal GW 
signal can be described by modelling the merged system as 
a perturbed Kerr SMBH [ringdown phase). The dynamics 
of the inspiral and ringdown phases are well understood and 
theoretical waveforms for the GWs have been derived. How¬ 
ever, certain parameters of the ringdown solution depend on 
unknown details of the merger phase and must be guided by 
the results of numerical simulation (see > 12 . 21 . 

Detection of GWs from a single binary source will be 
complicated by the presence of galactic and extra-galactic 
GW foregrounds, and by LISA’s sensitivity to GWs from all 
sky directions. Matched hlter template searches will there¬ 
fore be necessary to detect the GW signal from an individual 
event (Hughes et al. 2001). The SNR for a matched hlter de¬ 
tection (p) is dehned by the ratio of the coherently folded 
signal and noise powers: 
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where h{f) = e^'"^^^h[t)dt is the Fourier transform of 

the dimensionless strain and Sh{f) is the spectral power of 
the noise. We calculate SNRs for the hnal year (before the 
merger phase) of the inspiral signal. The effective duration 
of the ringdown phase is much shorter than one year and we 
therefore calculate ringdown SNRs over the entire phase. Es¬ 
timates of LISA’s noise spectrum are evolving with plans for 
the instrument design. A recently suggested sensitivity goal 
(Bender 2003) has a (dimensionless) threshold sensitivity 

dehned as , , for p = 5 and 1 yr of observation! 

^ V5x3.15x10’’s’ ^ ■' ■' 

with a low frequency power-law slope of —2.5 for frequencies 
ranging between 10~® and 10“^ Hz, and a power-law slope 
of —3 for frequencies of 3 x 10“® to 10“® Hz. Bender (2003) 
suggests a hard low frequency cut oh at /min = 3 x 10“® Hz, 
roughly corresponding to the lowest resolvable frequency for 
a year long mission. LISA will be sensitive to frequencies up 
to /max ~ 1 Hz^. As a general rule p 5 is required for the 
conhdent detection of a signal (Hughes et al. 2001). 


2.1 Inspiral phase 

As the orbital separation of the binary shrinks, both the am¬ 
plitude and frequency of the inspiralling strain increase. To 
calculate accurate inspiral phase SNRs, it is important to in¬ 
clude the frequency dependence of the strain in the integra¬ 
tion of signal to noise. We achieve this using the technique 
outlined by Flanagan & Hughes (1998, henceforth FH98). 


^ http://www.srl.caltech.edu/%7Eshane/sensitivity/ 
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FH98 show that sky averaged squared SNR ((p^)) can be 
re-expressed in terms of the emitted GW energy spectrum 

dEe 

dfe 
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where fe is the frequency of the GW. This is related to the 
observed GW frequency (/) via the cosmological redshift of 
the source [/ = /e/(l + z)]- The inspiral energy spectrum 
may be easily calculated to first order from the quadrupole 
approximation for the gravitational wave luminosity (FH98) 


equation 0 accounts for the reduced signal amplitude for 
unequal mass BBH coalescence (FH98). Numerical simula¬ 
tions suggest that ~ 0.01 — 0.03 (Baker 2001; FH98 and 
references therein). We adopt trd = 0.01. Assuming that 
the distribution of the strain between polarisations (h+, hx) 
mimics that of the inspiral phase, the dependence of h+ and 
hx on the orientation of the system’s angular momentum 
(L) and sky position [h{9, <j))] is given by (Fryer, Holtz & 
Hughes 2002) 

h+{t) = ^+exp(-7r/rd,obsl/Q)cos(27r/rd,obsf + A); 

hx{t) = ^xexp(-7r/rd,obsf/Q)sin(27r/rd,obsi + A), 


dEe _ 1_ 2/3^2/3 AfbhAMbh r-1/3 

dfe ~ s'" (Mbh + AMbh)i/3-^" ’ 

= (3) 

where Mhh and AMhh are the individual BBH masses corre¬ 
sponding to a reduced mass p = total mass 

Mtotai = Mbh + AMhh- Numerical simulations of BBH coa¬ 
lescence indicate that this approximation is reasonable un¬ 
til / reaches fmergeiz) = 205(1 -I-z)“^(Mtotai/20M©)“^ Hz 
(FH98) which roughly corresponds to the observed fre¬ 
quency when equal mass BBHs have an orbital separation 
of 3 Schwarzchild radii (). This frequency forms the 
upper limit of integration in equation 0. Since the GW 
frequency of BBH coalescence typically increases with time, 
this frequency may be used to delineate the inspiral and 
merger phases. The observed GW frequency of the inspiral 
signal at an observed time t' prior to the merger phase is 
given by (FH98) 


fit 
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Depending on the binary characteristics, /oneyr = f{t' = 
3.15 X lO^s) may lie below the LISA waveband. We there¬ 
fore use the greater of /oneyr and /min as the lower limit of 
integration in equation 0. 


where 

A+ = Ad[l -I- {L.hf]-, 

Ax ~ 2vTrd(L.?r), 

/rd,obs = 4he observed ringdown GW frequency and 

A is the initial phase. Leaver (1985), Echeverria (1989) and 
Fryer Holtz & Hughes (2002) have produced fitting formulae 
for /rd and Q in terms of Mtotai and a: 

/Ad = ^Hz(^^) [1-0.63(1(6) 

Q = (7) 

which are accurate to « 5 per cent. Hughes & Blandford 
(2003) note that the value of a will depend on the role 
that SMBH merging plays in SMBH evolution. We fol¬ 
low Hughes (2002) and adopt a = 0.997. The observed 
strain [H{t)] is the sum of the strain in each polarisation 
weighted by the corresponding detector response functions 
[F+{0, 4>, Ip), Ex {9, (p, Ip)', Thorne 1987] for a source with po¬ 
larisation axes rotated by an angle xp to h{9, <p) 


Hit) = E+i9, <p, xP)h+{t) + Ex (9, </, xP)hx it). (8) 

For a ~ 1, the quality factor is large and most of the 
energy radiated during the ringdown phase is observed at 
the central ringdown frequency (/rd.obs). In this limit the 
SNR may be approximated by (Hughes 2002) 


2.2 Ringdown phase 


We use the technique described by Hughes (2002) to cal¬ 
culate SNRs for the ringdown phase. BBH coalescence is 
expected to result in a rotating SMBH with a bar-like ex¬ 
citation of the event horizon (Hughes 2002). This distor¬ 
tion of the Kerr solution may be modelled by the bar-like 
[H, m) = (2,2)1 quasi-normal mode (Kokkotas & Schmidt 
1999). The corresponding time dependent strain is described 
by an exponentially damped sinusoid with total (observed) 
amplitude (.4rd) whose value depends on the fraction of the 
BBH mass-energy (Mtotaic^) radiated during the ringdown 
phase (crd), the binary mass, and the spin parameter of the 
merger product (a). The amplitude is 


Alrd = 


(^) 1 

/5erd Mtotai 

4p 1 

\cj Deiz)\ 
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where Deiz) is the comoving distance to a source at red- 
shift z. The quality factor (Q) is related to the damp-time 
of the GW signal (r) by Q = Tvrfrd- The final factor in 


2 2j°°H^it) 
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For comparison with the sky averaged inspiral SNR we av¬ 
erage the SNR in equation m over randomly generated sets 
of i9, (p, xp. A, L.h). 


2.3 Signal to noise ratios 

More massive BBH systems radiate GWs at lower frequen¬ 
cies and for a given binary system the ringdown signal is 
at a higher frequency [see equations @ & ®]- Therefore 
if the sensitivity curve has a low frequency cut off and/or 
degrades rapidly with frequency, searches for the ringdown 
signal will probe more massive BBHs. We find that ringdown 
searches can detect coalescence between BHs with masses 
up to Mbh ~ 10®’^ Mq (10®'^ M©) at 2 = 1 (z = 9). In 
the left hand panel of Figure Q we plot the SNR for the 
inspiral idot-dashed line) and ringdown isolid line) phases 
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Mbh (^^sol.r) 


Figure 1. Signal to noise ratios for the inspiral {dot-dashed line) and ringdown {solid line) phases of equal mass binary black hole 
coalescence as a function of single black hole mass (left panel). For completeness we show the ringdown signal to noise ratios for random 
sky angles and orientations {dots). On the right we plot the ratio of ringdown to inspiral signal to noise as a function of coalescing black 
hole mass. The hatched region denotes the mass range over which the ringdown signal allows confident detection {p ^ 5). 


of equal mass, 2 = 1 BBH coalescence as a function of 
single BH mass. The right hand panel shows the ratio of 
ringdown to inspiral SNRs, with the hatched region de¬ 
noting the mass range over which the ringdown signal is 
detectable. The SNR is greater in the ringdown signal for 
Mhh > 10® Mq . We note that there exists a range of SMBH 
masses (Mbh ~ 10®'^ — 10®'^ M© at 2 = 1) for which BBH 
coalescence is detectable only in the ringdown phase. Since 
characteristic frequencies decay with (I-I- 2 ), this mass range 
decreases with redshift. By 2 = 9, the mass range is compa¬ 
rable to the typical SMBH mass observed in the local uni¬ 
verse (Mbh ~ 10^'® — 10®'® Mq). These results are consistent 
with the calculations of FH98 and reinforce the point that 
searches for gravitational radiation from the ringdown phase 
may generate constraints on the high mass SMBH popu¬ 
lation. The higher frequency of the ringdown signal from 
SMBH binaries offers limited advantage for their detection 
if the low frequency sensitivity has a power-law slope of —2 
and no lower cut-off (corresponding to the optimal scenario 
of constant accelerometer dominated noise, see White 2002). 
However, it is expected that frequency dependent sources of 
spurious acceleration will result in a steeper low-frequency 
sensitivity curve (Bender 2003). 

In addition to the instrumental noise, unresolved GWs 
from various classes of binary systems are expected to limit 
LISA’s sensitivity (see Sesana et al. 2004b). However, Sesana 
et al. (2004b) calculate that BBH coalescence only domi¬ 
nates the sky noise below a frequency of Hz, and this 
noise is still an order of magnitude lower than the suggested 
LISA sensitivity. Therefore the unresolved BBH background 
should not significantly affect LISA’s sensitivity to super- 
massive BBH coalescence. 


3 MERGER AND COALESCENCE RATES 

Estimates of the merger rate of dark matter halos are pro¬ 
vided by the extended Press-Schechter (1974) formalism 
(Lacey & Cole 1993). We denote the merger rate per halo of 
mass M per unit time with halos of mass between AM and 
AM -I- dM by I m ■ The local (all-sky) detection rate 


of SMBH mergers (f^) occurring in the redshift interval 
(2,2 -b dz) may be estimated (e.g. Wyithe & Loeb 2003) by 

integrating the halo merger rate over the Press- 

Schechter dark matter halo mass function (J^) as below: 


d'N f°° A».r dn 

—— = dM dAM(e(Mbh,AMbh,2)):777 

dzdt Jq J m /3 dM 




dAMdt 


I M/3 

( 10 ) 


XT 

where is the comoving volume per unit solid angle be¬ 
tween 2 and z dz and the factor (1 -|- 2)“^ accounts for 
time-dilation. The product of the halo SMBH occupation 
fractions [/bh(M, 2)/bh(AM, 2)] relates the expected SMBH 
merger rate to the halo merger rate. 

Equation lHUl has the following additional features. We 
assume an observationally motivated dependence of SMBH 
mass on host halo mass [Mbh = Mbh(M, 2), see S 13.11 . The 
merger rate is weighted by the probability for detection 
of the central SMBH coalescence, denoted (0), where the 
detection 0 is determined through calculation of the SNR 


0 (Mbh, AMbh, 2 ) = 0,if p(Mbh, AMbh, 2 ) < 5, 

= l,if p(Mbh, AMbh,2) ^ 5. (11) 

We exclude halo mergers where the accreted satellite halo 
takes longer than the Hubble time to sink to the centre of 
the merged system, since the BBH will not have time to 
form. This is achieved by limiting the merger rate contribu¬ 
tion to halo pairs with mass ratios smaller than 3 (Colpi, 
Mayer & Governato 1999). We assume efficient coalescence 
of SMBH binaries in halo mergers, making our prediction 
of the coalescence rate a maximum. The merger rate scales 
with the fraction of BBHs that evolve into the GW domi¬ 
nated regime (emrg). We further assume coalescence to be 
rapid. This last assumption may affect the redshift distri¬ 
bution of our expected counts (although this effect may be 
small, see Sesana et al. 2004a), but should have little effect 
on estimates of the total event rate. 
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3.1 The Afbh—A^haio relationship 

Shields et al. (2003) investigated the evolution of the rela¬ 
tionship between Mhh and the velocity dispersion of the host 
bulge (cTc) out to z ~ 3 using a sample of radio-quiet AGN. 
They found no evidence for variation from the locally de¬ 
termined relationship (Ferrarese & Merritt 2000; Gebhardt 
2000) 


Mbh = 10®- 


200kms 


( 12 ) 


Ferrarese (2002) found that the velocity dispersion of the 
local galactic bulges correlates with the circular velocity in 
the fiat region of the halos rotation curve ctc ~ indicat¬ 
ing that Mbh ~ -Uc- Assuming that there is no evolution in 
the Oc-Vc relation, and relating Vc to the virial velocity of a 
virialised halo of mass M 


= 23.4 
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IOS/i-iMq, 


, 1/3 
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kms 


(13) 


[where flm (n^(i+z)ii+n^+n^( 1 + 2)2 , d _ 1 and Ac 

IStt^ + 82d — is the overdensity of a virialised halo at 
redshift z] we can determine the redshift dependent Mhh- 
Mhaio relationship (normalised to the local relationship from 
Ferrarese 2002) 


Mbh = lO^M© 


/ Mhaio 

5/3 

■l-b2:' 

V 1.5 X 1O12M0, 

) 
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(14) 


The Mbh“Mhaio relationship has been established lo¬ 
cally for halo masses in the range 10^^ — lO^^'^'M© (Fer¬ 
rarese 2002). Fits to gravitational lens separation distribu¬ 
tions [Porciani & Madau 2000; Kochanek & -White 2001 
hereafter KW2001] and simple models for baryon cool¬ 
ing (Gole et al. 2000; KW2001) indicate that halos with 
M > 3 X 10^^ M 0 tend to form groups and clusters of galax¬ 
ies. This will suppress the number of very massive SMBHs 
relative to the number of massive dark matter halos. In ad¬ 
dition, limitations on accretion from the IGM, and feedback 
from massive galactic winds are expected to suppress star- 
formation in small mass halos (see Barkana & Loeb 2001 and 
references therein), and may also inhibit SMBH formation. 
In the following section, we construct an observationally mo¬ 
tivated SMBH occupation fraction (/bh) that empirically ac¬ 
counts for these effects. 


3.2 The black hole occupation fraction 

A model which assumes that all dark-matter halos contain 
SMBHs (e.g. Wyithe & Loeb 2003) may substantially over¬ 
estimate the SMBH coalescence rate. Here we assume that 
all galaxies residing in halos with virial temperatures larger 
than a critical value host a central SMBH [with mass deter¬ 
mined by the Mbh-Mhaio relationship, equation Jill]. We 
construct the galaxy occupation fraction of dark-matter ha¬ 
los at z = 0 by comparing the observed velocity distribu¬ 
tion of local galaxies with a velocity function that is gen¬ 
erated from the Press-Schechter mass-function, and which 
assumes a relationship between cooled baryonic and virial 


velocities (see Motivated by the absence of evolu¬ 

tion in the observed correlation between Mbh and the veloc¬ 
ity dispersion of the host bulge, we assume that the fraction 
of halos of mass (M) and redshift ( 2 ) that host a central 
SMBH is determined by v^ir independently of redshift [ie. 

/bh(M, 2 ) = /bh(Wvir)]. 


3.^.1 Galaxy velocity functions and SMBH demography 

Measurements of the circular velocities of galaxies reflect the 
dynamics of the stars in the inner region of the halo where 
the density profile differs significantly from the profile pre¬ 
dicted by Navarro, Frenk & White (NFW, 1997) for GDM. 
Simple galaxy formation models consider the effect of adia¬ 
batic cooling of baryons on the radial density and velocity 
profiles of galaxies but neglect heating due to starformation. 
Adiabatic cooling and condensation of the halo’s baryonic 
component into a rotationally supported disk steepens the 
inner density profile, resulting in a higher maximal circular 
velocity (ucooi) than implied by the initial profile (Dalcanton, 
Spergel & Summers 1997; Mo, Mao & White 1998; Gonzalez 
et al. 2000; KW2001). 

To model the cooled halo density profile, Mo et 
al. (1998) assume that the baryonic and dark matter com¬ 
ponents are initially uniformly mixed in an NFW density 
profile with concentration parameter c, spin parameter A 
and specific angular momentum j = 0.05. The peak circular 
velocity of the modified density profile may be determined 
by self-consistently solving (with angular momentum and 
energy conservation) for the scale factor of the exponential 
disk. For a disk with specific angular momentum equal to 
that of the original halo, the cooled velocity is related to the 
virial velocity by (Mo et al. 1998) 

/ ^ \ -0.06+2.71md+0.0047/A 
Ucool « (1 - 3md -I- 5.2m^) f ^ j 

X (1 - 0.019c -b 0.0025c^-b0.52/c)uvir. (15) 

Numerical N-body simulations of hierarchical structure for¬ 
mation suggest the average halo spin to be roughly inde¬ 
pendent of halo mass and equal to about 0.04 (Bullock et 
al. 2001a). Furthermore, the concentration parameter of the 
NFW profile may be related to the halo virial mass Mvir (at 

2 = 0) by c = 9 (Bullock et al. 2001b). 

We assume the fraction of halo mass in the cooled disk to 
be rritj = 0.05, independent of halo mass. The cooled baryon 
fraction is thought to be smaller in more massive halos (see 
KW2001), however we find that allowing mj. to vary over 
halo mass as suggested by the cooling model of KW2001 
has little effect on the profile of the adiabatically cooled ve¬ 
locity function. The halo cooled velocity function may then 
be computed from the PSMF through the change of variable 

dn dn 

dVcaal dM 

For a cooled baryon fraction (md) that is independent of 
mass, this prescription results in a Mhaio ~ -^^cooi depen¬ 
dence of cooled velocity on halo mass for 10® M© < Mhaio < 
10^® Mq . Consequently, the net effect of this adiabatic cool¬ 
ing model is to translate (and re-scale) the Press-Schechter 


dM 


dVcool 


(16) 
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Figure 2. Occupation fraction [solid black line) of SMBHs in 
halos as a function of the maximum circular velocity of the cooled 
density profile (lower axis) or the halo virial velocity (upper axis). 
The solid grey and dashed grey lines show the (local) cooled Press- 
Schechter halo velocity distribution and observed galaxy velocity 
distribution (calibrated by the right hand axis) as described in 


distribution of halo virial velocities. To construct our oc¬ 
cupation fraction we assume that there is one circular ve¬ 
locity for which galaxies have an occupation fraction of 
unity. We then find the (linear) relationship between cooled 
and virial velocity that enforces this condition. We have 
adopted the observational galaxy velocity function used by 
KW2001 for this purpose. KW2001 provide data down to 
Wcooi ~ 70 km s“^. Below this we extrapolate with a sin¬ 
gle power law ~ with /3 = — 1 (appropriate for a 

luminosity function with faint-end slope —1). 

This process returns a halo cooled velocity (ucooi) which, 
at the peak of the occupation fraction, is consistent with the 
velocity derived using the prescription of Mo et al. (1998) to 
within « 15 per cent. Moreover, our relationship is con¬ 
sistent with the ratio of Uvir/ucooi ~ 1.8 measured from 
weak leasing studies (Seljak 2002). The resulting occupa¬ 
tion fraction (see Figure lU peaks at Ucooi ~ 230 km s ^ 
(uvir ~ 125 km s“^). This cooled velocity is also consis¬ 
tent with the circular velocity at which Press-Schechter and 
observed galaxy velocity functions agree in the works by 
Gonzalez et al. (2000) and KW2001. 

van den Bosch et al. (2002) provide another estimate of 
the dependence of galaxy formation efficiency on host halo 
mass via the average halo mass-to-light ratio. This quantity 
may be estimated by constructing halo mass conditional lu¬ 
minosity functions using the relationship between (luminos¬ 
ity dependent) observed galaxy-galaxy correlation lengths 
and mass-scale for a hierarchical cosmology (see van den 
Bosch et al. 2002). Their calculation yields a minimum mass- 
to-light ratio (indicative of maximum galaxy formation effi¬ 
ciency) for dark matter halos with Mhaio ~ 3 x 10^*^ M©, 
which is comparable to the peak in fhh{M,z = 0) at 
M ~ 4.8 X 10“h“i M©. 


'°giO C^bh.min / “solar) 

0 2 4 6 8 10 



Figure 3. Expected detectable binary black hole coalescence 
rates as a function of the minimum virial temperature of a black 
hole hosting halo (Tvir,mini lower axis) or the minimum black hole 
mass (Mbh.mini upper axis) for both inspiral {dashed line) and 
ringdown {solid line) phases assuming efficient coalescence. The 
dotted lines highlight the predicted event rates if the minimum 
SMBH mass is taken to be the minimum SMBH mass in local 
galactic bulges (M^h min ~ 10^ M©) or the typical mass of a 
SMBH powering a quasar at z ~ 2 (Mbh.min ~ 10^ M©). 


The occupation fraction decays rapidly for cooled veloc¬ 
ities larger than 230 km s“^. We note that errors in /bh at 
high Uvir will not effect our predicted event rates due to the 
exponential tail of the PSMF. Indeed Wyithe & Loeb (2003) 
found that imposing a sharp cut-off to /bh at a virial velocity 
of 600 km s“^ had little impact on their predictions. 


3.2.2 The minimum SMBH mass 

Prior to the reionisation of the IGM, halos required a virial 
temperature larger than ~ lO'* K for the accreted gas to cool 
via atomic hydrogen transitions within a Hubble time. After 
reionisation of the IGM, the higher IGM temperature inhib¬ 
ited the accretion of gas into halos with virial temperatures 
lower than ~ 2 x 10® K. These virial temperatures specify 
the minimum mass of a halo inside of which a SMBH may 
form (e.g. Barkana & Loeb 2001) 



Equivalently, a minimum virial temperature specifies a min¬ 
imum depth for the potential well (or velocity, or mass) of 
a halo in which a SMBH can exist. Given the Mbh-A7haio 
relationship, equation llTTt implies a minimum SMBH mass 
(Afbh.min) which is independent of redshift. For Tvir.min = 
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Figure 4. Redshift distribution of ringdown (left panel) and inspiral (right panel) phase detections of binary black hole coalescence 
for several Tvir^min assuming rapid, efficient coalescence. The total event rates (A'^total) for each Tvir min curve are listed in order of 
descending peak detection rate. 


10"^ K (2 X 10® K) this minimum SMBH mass corresponds 
to Mbh.min ~ 15 Mq (1.5 X 10^ Mq). In contrast, the lo¬ 
cally detected population of SMBHs reside in more massive 
halos with virial temperatures greater than « 10®'^ K, al¬ 
though this lower limit may be partially attributable to the 
difficulties in obtaining dynamical evidence for BHs with 
Mhh S 10^ (e.g. Kormendy & Richstone 1995). 

3.3 Event-rate predictions 

The predicted event rates are plotted in Figure |3 as a func¬ 
tion of the minimum virial temperature (lower axis) and 
minimum SMBH mass (upper axis). The event rates per 
year decline with increasing Tvir.min- For Tvir.min ~ 10® K 
(IHbh.min ~ 10"' M©) we predict « 40 detections per year 
in each of the ringdown and inspiral phases. These predic¬ 
tions are significantly below those of Wyithe & Loeb (2003), 
with the discrepancy due to our inclusion of the SMBH oc¬ 
cupation fraction in the merger rate. Our SMBH occupation 
fraction and accurate detection criteria result in predictions 
that decay more rapidly with increasing Tvir.min than previ¬ 
ous calculations suggested. 

SMBHs powering luminous quasars at 2 ~ 2 weigh 
> 10^ Mq. If Mbh.min ~ lO^®' Mq the predicted event- 
rates are very low; falling below 0.001 yr“' for each phase. 
Hence LISA will not detect coalescence from the major 
merges of bright quasar hosts. However, there is evidence for 
SMBHs as small as « 10® Mq in local galactic bulges (Fer- 
rarese 2002) and low-luminosity, low-redshift AGN (Greene 
& Ho 2004). A reasonable estimate for the minimum SMBH 
mass is therefore Mbh.min = 10® Mq, which corresponds to 
predictions of ~ 15 events per year in each phase. Figure 3] 
shows the redshift distribution of events as a function of the 
minimum virial temperature. Most of the events originate at 
z « 3 — 4 due to mergers between halos more massive than 
10^° Mq. 

Sesana et al. (2004a&b) make predictions for the de¬ 
tection rate of inspiral phase BBH coalescence for a seeded 
SMBH formation model; estimating « 35 “bursts” in 3 years 
of observation due to events at « « 4. This calculation 


employs a specific model for the formation and evolution 
of BBHs embedded in an isothermal stellar distribution. 
Sesana et al. (2004a) use the seeded SMBH growth model 
of Volonteri et al. (2003) in which (large) stellar mass BHs 
form in 3.5 — cr peaks at z = 20 and grow via accretion and 
coalescence during major halo mergers. 

While we also rely on the merger rate of dark matter 
halos for the basis of our coalescence rate, our calculation 
differs by employing empirically motivated estimates for the 
population of SMBHs in halos. We do not assume the physics 
leading to the absence of galaxies (and by extension SMBHs) 
in large and small halos, however we have shown that the 
inclusion of the occupation fraction significantly affects the 
coalescence rate. We have assumed an efficient coalescence 
of SMBH binaries, while Sesana et al. (2004a) specify a de¬ 
tailed model for binary evolution that results in some bina¬ 
ries being ejected from the galaxy core (where they do not 
coalesce). 

We estimate an upper limit for the rate of supermassive 
BBH coalescence (where both binary members weigh more 
than 10® Mq) of ~ 15 yr“', which is larger than, but compa¬ 
rable to, the ~ 9 yr“' estimated by Sesana et al. 2004a&b for 
binary coalescence with one member heavier than 10® Mq. 
If intermediate mass black holes exist, then our estimate for 
the detectable coalescence rate increases only modestly [due 
to the declining occupation fraction for lower mass halos (see 
Figure 121] to ~ tens per year (see Figure EJ. 


4 CONCLUSION 

We have computed the SNR of both the inspiral and ring- 
down phases of SMBH coalescence, and calculated mass 
and redshift ranges over which BBH coalescence will be de¬ 
tectable by LISA. The higher frequency of the ringdown sig¬ 
nal allows detection of the coalescence of more massive bi¬ 
nary systems. In particular, searches for gravitational waves 
for the ringdown signal offer a superior probe of BBH co¬ 
alescence between SMBHs which are more massive than 
« 10® Mq, and allow confident detection of the coalescence 
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of SMBHs in the mass range Mbh ~ 10®'^ — 10®’^ Mq at 
2 = 1 should they exist (or Mbh ~ 10^ ® — 10®’^ M© at 
2 = 9). 

We have predicted the rate of binary black hole coa¬ 
lescence detectable by LISA assuming a hierarchical cos¬ 
mology in which black holes form in halos with virial tem¬ 
perature above a critical value and obey an observationally 
motivated scaling with halo mass and redshift of the form 
Mbh ~ '^haio(^ + Prompted by the observation that 

the numbers of both large and small galaxies fall short of 
the expected numbers of dark-matter halos, we compute 
the empirical galaxy occupation fraction. This is obtained 
by comparing the locally observed galaxy velocity function 
with a cooled Press-Schechter (1974) velocity function. Our 
predicted event rate declines rapidly with increased esti¬ 
mates for the smallest SMBH that can form. If we are guided 
by the smallest nuclear BH mass inferred in local galactic 
bulges and nearby low-luminosity AGN (« 10® Mq), then 
we expect « 15 detections per year. Most of these events 
will be detectable via both their final-year inspiral and ring- 
down signals. Our observationally motivated estimates for 
the event rate of BBH coalescence therefore indicate that 
observations of coalescing BBHs are a realistic target for 
LISA. 
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